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Os ; Abstract 

a^ . 

Q ' The Fermi motion parameter Pp is the most important parameter of AC- 

D . CMM model, and the value Pp ~ 0.3 GeV has been used without clear theo- 

Q I retical or experimental evidence. So, we attempted to calculate the value for 

00 ' Pp theoretically in the relativistic quark model using quantum mechanical vari- 

CN . ational method. We obtained Pp ~ 0.5 GeV, which is somewhat larger than 

I , 0.3 GeV, and we also derived the eigenvalue of Eb — 5.5 GeV, which is in 

I reasonable agreement with nis = 5.28 GeV. We also recalculated |T4b/V"cfe| as 

' a function of p„. 

cn ; 

^ ■ The simplest model for the semileptonic i?-decay is the spectator model which 

^ i considers the decaying 6-quark in the 5-meson as a free particle. The spectator 

Q\ • model is usually used with the inclusion of perturbative QCD radiative corrections. 

^ I The decay width of the process B Xqlu is given by 

^' Tb{B XJu) ~ T,{b qlu) = \V,f{^^)f{'^)[l - --^(^)], (1) 

! where rriq is the mass of the g-quark decayed from 6-quark. The decay width of the 



spectator model depends on mf, therefore small difference of mi, would change the 
decay width significantly. 

Altarelli et al. ^ proposed their ACCMM model for the inclusive i?-meson semilep- 
tonic decays. This model incorporates the bound state effect by treating the 6-quark 
as a vitual state particle, thus giving momentum dependence to the 6-quark mass. 
The virtual state 6-quark mass W is given by 

W\p) = m| + m% - 2mB^V^ + m% (2) 

in the i?-meson rest frame, where ntsp is the spectator quark mass, uib the 5-meson 
mass, and p is the momentum of the 6-quark. 

^Talk given by D.S. Hwang, at the Int. Workshop on B-physics, Nagoya, Japan, on Oct. 26-28, 
1994. Proceedings pubhshed by World Scientific, Singapore, edited by A. Sanda and S. Suzuki 
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For the momentum distribution of the virtual 6-quark, Altarelh et al. considered 
the Fermi motion inside the 5-meson with the Gaussian momentum distribution 



e-P'/P^, (3) 



where the Fermi motion parameter is treated as a free parameter. And the decay 
width is given by integrating the width Ffe in (1) with the weight 0(p). Then the 
lepton energy spectrum of the 5-meson semileptonic decay is given by 

-T^{Pp,msp,mq,mB) = dp p (j){p) —{mb = W,mg}, (4) 
(Wl J tt-C/i 

where Pmax is the maximum kinematically allowed value of p = |p|. The ACCMM 
model, therefore, introduces a new parameter p^ for the Gaussian momentum distri- 
bution of the 6-quark inside 5-meson instead of the 6-quark mass of the spectator 
model. In this way the ACCMM model incorporates the bound state effects and 
reduces the strong dependence on 6-quark mass in the decay width of the spectator 
model. 

The Fermi motion parameter p^ is the most essential parameter of the ACCMM 
model as we see in the above. However, the experimental determinations of its value 
from the lepton energy spectrum have been very ambiguous until now because various 
parameters, such as p^, rrig and m^p, are fitted all together from the lepton energy 
spectrum, and because the perturbative QCD corrections are sensitive in the end point 
region of the spectrum. We think that the value p^ ~ 0.3, which has been widely used 
in experimental analyses, has no theoretical or experimental clean justification, even 
though there has been recently an assertion that the BSUV model [|| is approximately 
equal to ACCMM model at p^ ~ 0.3. Therefore, it is strongly required to determine 
the value of p^ more firmly when we think of the importance of its role in experimental 
analyses. The better determination of Pp is also interesting theoretically since it has 
its own physical correspondence related to the Fermi motion inside i?-meson. In this 
context we are going to determine theoretically the value of Pp in the relativistic 
quark model using quantum mechanical variational method. 

We consider the Gaussian probability distribution function 0(p) in (3) as the 
absolute square of the momentum space wave function x(p) of the bound state B- 



m = 47r|x(p)r, X(P) = y^jsT^^"' (5) 

The Fourier transform of x(p) gives the coordinate space wave function tp{r), which 
is also Gaussian, 

\/tt 



Then we can approach the determination of Pp in the framework of quantum me- 
chanics. For the i?-meson system we treat the 6-quark non-relativistically, but the u- 
or rf-quark relativistically with the Hamiltonian 

H = M + ^ + ^p^ + m^ + V{r), (7) 
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where M = rrif, is the 6-quark mass and m = rrisp is the u- or d- quark mass. We 
apply the variational method the Hamiltonian with the trial wave function 

V(r) = (^)3/2e-M^-V2^ (8) 



where n is the variational parameter. The ground state is given by minimizing the 
expectation value of H, 

{H) = {^p\H\^lJ) = Eifx), = at fx = ft, (9) 

and then p. = Pp and E = E{ji) approximates uib- The value of /i or corresponds to 
the measure of the radius of the two body bound state as can be seen from = 

and (r^)^ = |-. 

In we take the Cornell potential which is composed of the Coulomb and linear 
potentials, 

V{r) = + Kr. (10) 

For the values of the parameters etc (= fos), and the 6-quark mass rrib, we use 
the values given by Hagiwara et al. 

ac = 0.47 {as = 0.35), K = 0.19 GeV^, nib = 4.75 GeV, (11) 

which have been determined by the best fit of the (cc) and (bb) bound states. For 
comparison we will also consider ac = 0.32 (a^ = 0.24), which corresponds to as{Q^ = 
ml). 

Before applying our variational method with the Gaussian trial wave function 
to the S- meson system, let us check the method by considering the T{bb) system. 

The Hamiltonian of the T{bb) system can be approximated by the non-relativistic 
Hamiltonian 

H ^2mb + — + V(r). (12) 
rrib 

With the parameters in (pT]) (or with ac = 0.32), our variational method with the 
Gaussian trial wave function (||) gives Pp = fi = 1.1 GeV and E = E{fi) = 9.49 GeV. 
Here Pp = 1.1 GeV corresponds to the radius R{T) = 0.2 fm, and E{T) = 9.49 GeV 
is within 0.3 % error compared with the experimental value E^xp = = 9.46 GeV. 
Therefore, the variational method with the non-relativistic Hamiltonian (0) gives 
fairly accurate results for the T ground state. 

However, since the u- or d- quark in the i?-meson is very light, the non-relativistic 
description can not be applied to the i?-meson system. For example, when we apply 
the variational method with the non-relativistic Hamiltonian to the 5-meson, we get 
the results 

Pp = 0.29 GeV, E = 5.92 GeV for a, = 0.35, (13) 
Pp = 0.29 GeV, E = 5.97 GeV for a, = 0.24. (14) 



The above masses E are much larger compared to the experimental value rriB = 5.28 
GeV, and moreover the expectation values of the higher terms in the non-relativistic 
perturbative expansion are bigger than those of the lower terms. Therefore, we can 
not apply the variational method with the non-relativistic Hamiltonian to the B- 
meson system. 

Let us come back to our Hamiltonian (0) of the 5-meson system. In our variational 
method the trial wave function is Gaussian both in the coordinate space and in 
the momentum space, so the expectation value of H can be calculated in either 
space from {H) = {ip{r)\H\ip{r)) = {x{p)\H\x{p))- Also, the Gaussian function 
is a smooth function and its derivative of any order is square integrable, thus any 
power of the Laplacian operator is a hermitian operator at least under Gaussian 
functions. Therefore, analyzing the Hamiltonian (|^) with the variational method 
can be considered as reasonable even though solving the eigenvalue equation of the 
differential operator (|^) may be confronted with the mathematical difficulties because 
of the square root operator in (|^). 

With the Gaussian trial wave function (j^) or (^, the expectation value of the 
Hamiltonian (|^) can be calculated easily besides the square root operator, 

(p^) = Mr)|p2|^(r)) = (x(p)|p^|x(p)) = ^/i^ (15) 
{V{r)) = (V;(r)|-^ + irr |^(r)) = — (-a,/i + ir/^). (16) 



Now let us consider the expectation value of the square root operator in the momen- 
tum space 



^ Jo 



oo 



■"Yx2 + (m//i)2 x'rfx. (17) 



The integral ( [T7| ) can be given as a series expansion by the following procedure. First, 
define 

/■oo rco 

lis) ^ / \G;^sx^e-'' dx = s'^ Vt^TT t V^* rft, (18) 
Jo Jo 

/* oo /* OO 

Jo(s) = / V^^!^T^e-^'dx = s VWTle-'''dt. (19) 
Jo Jo 

Next, from (|l^) and ([l9|) , we find the following differential relations 

Combining two equations in (pOD, we get a second order differential equation for /(s), 

sl"{s) - (1 + .)/'(.) + \l{s) = 0. (21) 



The series solution to (|2T|) is given as 

lis) 



his) 
his) 



ci/i(s) + 02/2(5), 
s^F{s;^,3) = s 



s H H 

2 32 192 



s' + 



1024 



5' + ---}, (22) 



se 



[his)? 
1/1 5 
- ^(1 + 2^+32 



(is 



16 



-s^lns 1 



2 7 

H 

32 192 



s' + 



1536 



-s H 

2 32 



s^ + - 



where -F(s; |, 3) is the confluent hypergeometric function which is convergent for any 



finite s, and the integral constants Ci ~ —0.095, C2 
derivation of these numerical values for c, . 



■1. See Appendix for the 

Finally, collecting (p^, (|16]) and ([T7|), the expectation value of H is written as 
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M + — fV 



H ^(-ftcAt + -^/yw) 



TT 



2/i 



(32 ~ 2ci)(m//i)^ + -(m//i)^ln(m/^) 



(23) 



up to (m/ /i)^ 



With the input value of m 
we obtain 



m 



sp 



0.15 GeV, we minimize (if) of (PB|), and then 



0.54 Ge\/, 
0.49 GeV, 



TUB 



-- E 
E -- 



■ 5.54 GeV 
5.63 GeV 



for 

for a. 



0.35, 
: 0.24. 



(24) 



Here let us check how much sensitive our calculation of Pp is by considering the case 
for comparison. For m^p = the integral in (p!7| ) is done easily 



where m = iris. 

and we obtain the following values of fl = by the above variational method. 



0.53 GeV, 
0.48 GeV 



E 
E 



5.52 GeV 
5.60 GeV 



for as 
for a. 



0.35, 
0.24. 



(25) 



As we see in (25), the results are similar to those in (p^ where m^p = 0.15 GeV. We 
could expect this insensitivity of the value of p^ to that of m^p because the value of 
rrisp, which should be small in any case, can not affect the integral in (23) significantly. 



The calculated values of the i?-meson mass, E, are much larger than the measured 
value of 5.28 GeV. The large values for the mass are originated partly because the 
Hamiltonian (psp does not take care of the correct spin dependences for B and B*. 
The difference Detween the pseudoscalar meson and the vector meson is given arise 
to by the chromomagnetic hyperfine splitting, which is given by 



V 



a. 



(26) 



Then the expectation values of Vg are given by 

2 ttc/i^ 



(Vs) 



for B, (K) 



for B' 



^ Mm ~" ' Mm 
and we treat (Vg) only as a perturbation. Then we get for B meson 



and for B* 



Pf 
Pf 



Pf 
Pf 



0.54 GeV, 
0.49 GeV, 



0.54 GeV, 
0.49 GeV, 



Eb = 5.42 GeV 
Eb = 5.56 GeV 



Eb* 
Eb* 



5.58 Ge\/ 
5.65 GeV 



for = 0.35, 
for as = 0.24, 



for as 
for a. 



0.35, 
0.24. 



(27) 



(2^ 



(29) 



The calculated values of the 5-meson mass, 5.42 GeV (a^ = 0.35) and 5.56 GeV 
{as = 0.24) are in reasonable agreement compared to the experimental value of = 
5.28 GeV; the relative errors are 2.7% and 5.3%, respectively. However, for Fermi 
motion parameter p^, the calculated values, 0.54 GeV (a^ = 0.35) and 0.49 GeV 
{as = 0.24), are somewhat larger than the value 0.3 GeV, which has been widely 
used in the experimental analyses of energy spectrum of semileptonic i?-meson decay. 
The value Pp = 0.3 GeV corresponds to the i?-meson radius Rb ~ 0.66 fm, which 
seems too large. On the other hand, the value p^ = 0.5 GeV corresponds to Rb ~ 0.39 
fm, which looks in reasonable range. 

If we use Pp = 0.5 GeV, instead of p^ = 0.3 GeV, in the experimental anal- 
ysis of the end point region of lepton energy spectrum, the value of iKifc/Kbl be- 
comes significantly changed. Using the experimental data of the end point region, 
z.e. 2.3 GeV < Ei < 2.6 GeV of the CLEO resuh M, we can find the relation 



Vub 


2 


Vub 


Vb 


p„=0.5 


Vcb 



=0.3 



£(03) 



V 



ub 



V 



cb 



X 1.81, 



=0.3 



(30) 



where T{pp) = dEi-^^{pp) with \Vub\ = 1, and we assume the value of \Vcb\ 
is determined independently from the other analyses. We numerically calculated 
f(0.3)/f(0.5_) by using §) with m,p = 0.15 GeV, rriq = 0.15 GeV and tub = 5.28 
GeV, to get its value as 1.81. Previously the CLEO analyzed with p^ = 0.3 GeV the 
end point lepton energy spectrum to get [ffl 



10' X IKfe/KbP 



0.57 ±0.11 (ACCMM [1]) 
1.02 ±0.20 (ISGW §). 



(31) 



As can be seen, those values are in large disagreement. However, if we use Pp = 0.5 
GeV, the result of the ACCMM model becomes 1.03 from (pO]), and these two models 
are in good agreement for the value of \Vub/Vcb\- Finally we show the values of 
\Vub{pF)/^ub{PF = 0.3)1 as a function of Pp in Fig. 1. 
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Appendix 

The integration constants Ci and C2 in (p2D are given by the following relations. 



1 r°° 1 
/(O) = --C2= / x'e-^ dx = -, (32) 
Z Jo 2 

I"(s^O) = 2ci + C2(-ilns-— ) 

= — / x\x^ + s)~^/^e-^ dx at s^O. (33) 
4 JO 

Then, from (0), we get 

C2 = -1. (34) 
The integral in (^) can be expanded as 

Jis = a^) = r x\x^ + a'')-^/''e-^"dx 
Jo 

x^[(x + ay - 2ax\~^''^e~^'' dx 

-3/2 







00 

x^ix + a)~^ 



2ax 



e ""^dx 



{x + ay 

^ (2n + l)!a" 

^0 2«(n!)2 Jo (x + a)2«+3^ y'^^) 
Next the integral in (|H3|) is obtained by successive differentiations of an integral. 



-e rfx = — — - — / e rfx. (36) 



(x + a)2"+3'' (2n + 2)! l^Sayl 



X + a 



Again the integral in ( ^6]) is related to another integral, for a small value of a, 

e - = 5] _^ ! + (_a)"+2 / -—dx. (37) 



X + a ^ 2V2/ Jox + a 



The integral in (|3^) can be expanded in a similar way as the series (|2^) was obtained 



by making use of differential equations. For a small value of a, 

/ _^dx = {2\na+^W + -- + ---) + V^e-^ a+7Ta' + 77^ + - ■ ■ , 38 

Jo X + a 2 2 21 3 5 2! 

where 7 ~ 0.5772 is the Euler's constant. In this way the constant ci is given by an 
infinite series, 

371^1 

ci = + — - - y ^ -0.0975. 39 

64 16 8 f"in2" ^ ' 
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Fig.l \VubiPp)/VubiPp = 0.3)1 as a function of p^,. 
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